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Abstract 

In the standard model the two-body charmless hadronic meson decays can occur via 
annihilation diagrams only. In this work, we studied the B^ PP, PV/VP, VV decays by 
employing the perturbative QCD (pQCD) factorization approach. From our calculations, we 
find that (a) the pQCD predictions for the branching ratios of the considered Be decays are in 
the range of 10"^ to 10"^ (b) for Be PV/VP,VV decays, the branching ratios of AS" = 
decays are much larger than those of AS" = 1 ones because the difi^erent Cabibbo-Kobayashi- 
Maskawa(CKM) factors are involved; (c) analogous to -B ^ Krj^'^ decays, we find Br(Bc 
K^rj') ~ 10 X Br{Bc K^rj), which can be understood by the destructive and constructive 
interference between the r]q and rjg contribution to the Be — > K^rj and Be — > K^rj' decay; (d) 
the longitudinal polarization fractions of Be VV decays are in the range of 86% — 95% and 
play the dominant role; and (e) there is no CP-violating asymmetries for the considered Be 
decays because only one type tree operators involved. 
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I. INTRODUCTION 



In 1998, a new stage of physics began because of the first observation of the meson 
Be at Tevatron [lij. For Be meson, one can study the two heavy flavors b and c in a 
meson simultaneously. From an experimental point of view, more detailed information 
about Be meson can be obtained at the Large Hadron Collider (LHC) experiment. The 
LHC is scheduled to start to run in this month, where the Be meson could be produced 
abundantly. The Be meson decays may provide windows for testing the predictions of the 
standard model(SM) and can shed light on new physics(NP) scenarios beyond the SM. 

From a theoretical point of view the non-leptonic decays of Be meson are the most 
complicated decays due to its heavy-heavy nature and the participation of strong inter- 
action, which complicate the extraction of parameters in SM, but they also provide great 
opportunities to study the perturbative and nonperturbative QCD, final state interac- 



tions, 
Refs. 
27, 



etc. The non-leptonic Be weak decays have been widely studied for example in 




factorization approach(QCDF) [34|, the perturbative QCD (pQCD) approach jss], 36, 13 
and other approaches and/or methods. 

In this paper we focus on the two-body non-leptonic charmless decays Be — >■ 
PP, PV/VP,VV (here P and V stands for the light pseudo-scalar and vector mesons), 
which can occur through the weak annihilation diagrams only. The size of annihilation 
contributions is an important issue in B physics. Indeed, the two-body charmless Be 
decays considered here are rather different from those Be J/ipPiV) decays where the 
initial c quark behaves spectator. 

Recently, the two-body non-leptonic charmless Be — > M2M3 ^ decays have been studied 
by using the SU (3) flavor symmetry or by employing the QCD factorization approach (38| . 
The authors in Ref. [3^ provided two different estimates for non-leptonic charmless Be 

decays. But their predictions for the branching ratios of Be — * <t>K^ ■, 'K* decays in the 
QCDF are much smaller (a factor of 10) than those obtained by using the SU(3) flavor 
symmetry. So large discrepancies among the theoretical predictions for the branching 
ratios indicate clearly that it is very necessary to make more studies for these kinds of 
Be decays by employing other different approaches, in order to understand these decays 
better and provide the theoretical support for the related experimental studies. 

In this paper, we will calculate the branching ratios and the polarization fractions 
of thirty Be — > PP, PV/VP, VV decays by employing the low energy effective Hamilto- 
nian [39] and the pQCD factorization approach. By keeping the transverse momentum 
kT of the quarks, the pQCD approach is free of endpoint singularity and the Sudakov 
formalism makes it more self-consistent. It is worth of mentioning that one can do the 
quantitative calculations of the annihilation type diagrams in the pQCD approach. The 
importance of annihilation contributions has already been tested in the previous predic- 
tions of branching ratios of pure annihilation B — >■ DsK decays [13], direct CP asym- 



metries of P° — s> vr+TT, K^-K decays |35|, |36|, |4l| and in the explanation of P — > (pK 



polarization problem 4^, 43], which indicate that the pQCD approach is a reliable method 



^ For the sake of simplicity, we will use M2 and M3 to denote the two final state light mesons respectively, 
unless otherwise stated. 
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to deal with annihilation diagrams. 

The paper is organized as follows. In Sec. [TTl we present the formalism and wave func- 
tions of the considered Be meson decays. Then we perform the perturbative calculations 
for considered decay channels with pQCD approach in Sec. IIIII The numerical results 
and phenomenological analysis are given in Sec. IIV[ Finally, Sec. |V] contains the main 
conclusions and a short summary. 



II. FORMALISM AND WAVE FUNCTIONS 



A. Formalism 



Since the b quark is rather heavy, we work in the frame with the Be meson at rest, i.e., 
with the Be meson momentum Pi = -^(1, 1,0t) in the hght-cone coordinates. For the 
non-leptonic charmless B^ — ^ M2M3 decays, assuming that the M2 (M3) meson moves in 
the plus (minus) z direction carrying the momentum P2 (P3) and the polarization vector 
€2 (e3)(if M2{3) are the vector mesons). Then the two final state meson momenta can be 
written as 

P2 = ^{l-rlrl0r), P, = ^{rll - rlOr), (1) 

respectively, where r2 = mM2/''^B,and r^, = niMz/f^B- When M2,M3 are the vector 
mesons, the longitudinal polarization vectors, €3 and eg , can be given by 

^ ■(l-r3^-r2^0^), eli = ^^{-rll-rlM. (2) 



The transverse ones are parameterized as = (0, 0, It), and ef" = (0, 0, It). Putting the 
(hght-) quark momenta in Be, M2 and M3 mesons as ki, /c2, and ks, respectively, we can 
choose 

ki = (XiP+,0,kiT), k2 = (x2P2^,0,k2T), k^ = {0, x^P-^ Mt) ■ (3) 

Then, for Be M2M3 decays, the integration over k^ , fc^, and kf will conceptually lead 
to the decay amplitudes in the pQCD approach, 

A{Be ^ M2Ms) ~ Jdxidx2dx3bidbib2db2b3db3 

•Tr [C{t)^BSxu bi)^M, {X2, &2)$M3 {X3, bs)H{xi, h, t)St{xi) e-^W|4) 

where 6, is the conjugate space coordinate of kix, and t is the largest energy scale in 
function H{xi, bi, t). The large logarithms ln{mw /t) are included in the Wilson coefficients 
C{t). The large double logarithms (In^ Xi) are summed by the threshold resummation 44 1. 
and they lead to St{xi) which smears the end-point singularities on Xj. The last term, 
Q-S(,t) ^ is the Sudakov form factor which suppresses the soft dynamics effectively [45[. Thus 
it makes the perturbative calculation of the hard part H applicable at intermediate scale, 
i.e., rriB^ scale. We will calculate analytically the function H{xi,bi,t) for the considered 
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decays at leading order(LO) in expansion and give the convoluted amplitudes in next 
section. 

For these considered decays, the related weak effective Hamiltonian i^e// iloj can be 
written as 

H^ff = ^ [V:,V^n {C,{fi)0,{f,) + CMO.if,))] , (5) 
with the single tree operators, 

O2 = Upl^il - 75)/^/3Ca7''(l - 75)&a , (6) 

where Vcb, Vud are the CKM matrix elements, "D" denotes the light down quark d or s and 
Ci{fi) are Wilson coefficients at the renormalization scale fi. For the Wilson coefficients 
Ci,2(At), we will also use the leading order (LO) expressions, although the next-to- leading 
order calculations already exist in the literature [39(]. This is the consistent way to cancel 
the explicit n dependence in the theoretical formulae. For the renormalization group 
evolution of the Wilson coefficients from higher scale to lower scale, we use the formulae 
as given in Ref. [36| directly. 



B. Wave Functions 



In order to calculate the decay amplitude, we should choose the proper wave functions 
of the heavy and light mesons. In principle there are two Lorentz structures in the -Bu,d,s 
or Be meson wave function. One should consider both of them in calculations. However, 
since the contribution induced by one Lorentz structure is numerically small 46|, |47| and 



can be neglected approximately, we only consider the contribution from the first Lorentz 
structure. 



(7) 



Since B^ meson consists of two heavy quarks and m^^ ~ 1^1, + rric, the distribution 
amplitude (pBc would be close to 6{x — rric/mBj in the non-relativistic limit. We therefore 



adopt the non-relativistic approximation form of as |19L l28 



■.6{x - mc/rriB,) 



where /^^ and are the decay constant of B^ meson and the color number, respectively. 
For the pseudoscalar meson(P), the wave function can generally be defined as. 



$p(x) 



:75{f0^(x) +m^0^(x) +m, 



lWp{x)] 



al3 



(9) 



where 



lA,P,T 



and rriQ are the distribution amplitudes and chiral scale parameter of the 
pseudoscalar mesons respectively, while x denotes the momentum fraction carried by 
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(a) (b) (0) (d) 



FIG. 1: Typical Feynman diagrams for two-body non-leptonic charmless Be decays. 

quark in the meson, and n = (1, 0, 0^) and v = (0, 1, 0^) are dimensionless light-like unit 
vectors. 

For the wave functions of vector mesons, one longitudinal(L) and two transverse(T) 
polarizations are involved, and can be written as, 

^vi^) = ^jMvrv'^M^) + rvMi^) + Mv<j)U^)}^^ , (10) 
^lix) = -^{Myi1^ct^l{x)+ilJrfv{^)+Mvit,.p.l,l'e*.^ ,(11) 

where ey denotes the longitudinal (transverse) polarization vector of vector mesons, 
satisfying P ■ e = in each polarization. We here adopt the convention e°^^^ = 1 for 
the Levi-Civita tensor e^'^"'^ . For the distribution amplitudes of pseudoscalar (f)p ' , 
and longitudinal and transverse polarization, (j)y^ and (py which will be presented in 
Appendix \M 

III. PERTURBATIVE CALCULATIONS IN PQCD 

From the effective Hamiltonian ([5]), there are 4 types of diagrams contributing to 
the Be M2M3 decays as illustrated in Fig. [T|, which result in the Feynman decay 
amplitudes pj^^^ and M^^^^, where the subscripts fa and na are the abbreviations 
of factorizable and non-factorizable annihilation contributions, respectively. Operators 
O12 are {V — A){V — A) currents, we therefore can combine all contributions from these 
diagrams and obtain the total decay amplitude as, 

A{Be ^ M2M3) = {fBF^l'''''a, + M^J^'^C,} , (12) 

where ai = Ci/3 + C2- In the next three subsections we will give the explicit expressions 
of Fj^^Ma^ M^^Ma ^^le decay amplitude A^B^ M2M3) for ^ M2M3 decays: 
including eight Be —>■ PP, fifteen B^ PV or B^ — > VP, and seven B^ — > VV decay 
modes. 

A. PP decays 

In this section, we will present the factorization formulae for eight non-leptonic charm- 
less Be —>■ PP decays. From the first two diagrams of Fig. [H i.e., (a) and (b), by per- 
turbative QCD calculations, we obtain the decay amplitude for factorizable annihilation 
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contributions as follows 







X {^/a(l - 2^3, 2^2, ^3, h2)Efa{ta) [x20^(x2)</)3 (Xg) + 2rlrl(p^{x^) 
X {{X2 + (0:2) + (2:2 - 1)02 (3^2))] + hfa{x2, 1 - X3, 62, h)Efa{tb) 

X [(X3 - 1)0^(X2)0^(X3) + 2roVoVf (0:2) ((X3 - 2)03^ (X3) - X30r(x3))] } ,(13) 

2f 3) 

where 02(3) corresponding to the distribution amplitudes of mesons M2(3), = 
^M2(Af3)^^^^^ and Cf = 4/3 is a color factor. In Eq. f|T3l) . the terms proportional to 
(rQ*-^-*)^ have been neglected because they are small indeed, max(rQ*-^'*)^ < 7%. The func- 
tion hfa, the scales and Ejait) can be found in Appendix [Bl 

For the non-factorizable diagrams (c) and (d), all three meson wave functions are 
involved. The integration of 63 can be performed using S function 5(63 — 62), leaving only 
integration of 61 and 62- The corresponding decay amplitude is 

M^^ = '^nCFTn'^Bc [ dx2dxs f bidbib2db2 

3 "Jo Jo 

X {/i^„(x2, X3, bl, b2)Ena{tc) [{tc - + l)<p^ {x2)<j)i {x^) + r^rjj (0f (0:2) 
X ((3rc + X2 - X3 + l)0f (X3) - (rc - X2 - X3 + l)0f (X3)) + 0^(x2) 

X ((rc - X2 - X3 + l)0f (X3) + (rc - X2 + X3 - 1)03 (0:3)))] - Ena{td) 
X [{n + re + X2 - 1)0^(X2)03 (X3) + ^o'^O (0f (2^2) ((4r6 + re + X2 - X3 
-l)0f (X3) - (re + X2 + X3 - 1)0![(X3)) + 0^(x2)((re + X2 + 0:3 - 1) 
X0f(a;3) - (re + X2 -X3 - 1)03(X3)))] /lKx2,X3,&l,&2)} , (14) 

where rf,(e) = mb(c)/mB,. 

For the 1] — 77'" system, there exist two popular mixing basis: the octet-singlet basis 



and the quark-flavor basis [48|, |49j. Here we use the quark-flavor basis 48|] and define 

rjg = [uu + dd)/\/2, rjs = ss. (15) 
The physical states rj and 77' are related to riq and ris through a single mixing angle 0, 

V \ _ TT(A\ fvA - f -sin0\ / riq 



?7 y \ils J y sm cos (p J \ T]s J 

We assume that the distribution amplitudes of 77^ and rjs are the same as the distribution 
amplitudes of vr, except for the different decay constants and the chiral scale parameters. 
The three input parameters fq, f^t and in the quark-flavor basis have been extracted 



from various related experiments [48|, |4^ 

/, = (1.07 ± 0.02)/^, /, = (1.34 ± 0.06)A, = 39.3° ± 1.0°. (17) 
The chiral enhancement factors are chosen as 



,2 



m 1 2 2 i , 2 • 2 i ^J^fs f 2 2\ I ■ 

rriQ = - — = - — m„ cos + m„, sm — [m, - mj cos sm ( 

2mg 2mq ' Jq 

^2 



mn" = ::r^ = [m^, cos^ + rri^ sin^ (m^/ — m^) cos0sin0]. (19) 



2m, 2m/ " ^ ^ ^2/, 



In the numerical calculations, we will use these mixing parameters as inputs. It is worth of 
mentioning that the effects of possible gluonic component of r( meson will not considered 
here since it is small in size 50, 51, 5^ - 

Based on Eqs. (fT2!) and fll3fl4l) . we can write down the total decay amplitudes for eight 
Be — > PP decays easily. 



TT TT 



-[/Sei^/a' Ci]|=0 

+ F;r^ai + M„^|^-^Ci]} COS0 
+ \iBFff\^ + Mir^C^\ I sin 



(20) 

(21) 

(22) 
(23) 



i^+TrO) 



A{B, K\+) 
A{B, ^ ir+77) 



^/a COS. 



F;^^%in, 



+ 



M„^>cos0-M„^j 



smi 



A{B, ^ ir+r;') = /iJ, 



F;>sin0 + F;, 



COS( 



+ 



Mt^^^ sin + Mlif^ cos 0] Ci } 



ai 



ai 



(24) 
(25) 

(26) 
(27) 
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B. PV, VP decays 

By following the same procedure as stated in the above subsection, we can obtain the 
analytic decay amphtudes for Be — > PV, VP decays, 

/■I roo 

" Jo Jo 

X {hfa{l - X3, X2, 63, b2)Efa{ta) [^^sc/'s (a;2) 03 (^^s) " 2ror3(/)^ (X3) 

X {{x2 + (a;2) + {x2 - 1)02 (2^2))] + hfa{x2, 1 - 2:3, 62, h)Efa{tb) 

X [(,T3 - l)0^(x2)03(a;3) - 2rlr^^P{x2) {{x^ - 2)0^(a;3) - a;30*(a;3))] } ,(28) 

-^tTo^ = ^^y^ 7rCirm|^ / dx2dxz [ bidbib2db2 

X {/i;^„(a;2, 0:3, 61, b2)Ena{tc) [{Tc -X2, + 1)4>2{X2)(I)Z{X^) - tIt^ {4h {^2) 

x((3rc + X2-X:i + 1)03(^^3) - (^c - a;2 - X3 + 1)03(2^3)) + 0^(2:2) 
x((rc - X2 - X3 + l)03(a;3) + (r^ - ^2 + ^3 - 1)0* (2:3)))] - Ena{td) 

X [(n + rc + X2- 1)02 (X2)03(.«3) ^ ^0^3 (0^ (a;2) ( (4rb + Tc + X2 - X3 

-l)0^(a;3) - (r, + 0:2 + 0:3 - l)03(x3)) + 0[(a;2)((rc + X2 + X3 - 1) 
x0^(a;3)-(r, + X2-X3-l)03(a;3)))]<a(a;2, 0:3,61,62)} , (29) 

pl poo 

FJJ[' = STrCpm^^ / dx2dxz / b2db2b^db2, 
" Jo Jo 

X {/i/a(l - 2:3, a;2, 63, b2)Efa{ta) [a;202 (2:2) 0^(2:3) + 2r2ro0f (a;3) 

x {{x2 + 1)02(2:2) + (a;2 - 1)02(0:2))] + hfa{x2, 1 - X^, 62, b^)Efa{tb) 

X [(;2:3 - 1)02(X2)0^(X3) + 2r2ro^0^(a:2) ((0:3 - 2)0f (X3) - x^cfi^ix;^ ] ,(30) 

-H^^ = — ^TrC^m^^ / dx2dx3 ( 61^6162(^62 

X {^na(3;2, 2:3, 61, b2)Ena{tc) [{Tc - + 1)02 (^2)0^ (Xg) + r2r|j (03(2:2) 

x((3rc + 2:2 - 2:3 + l)0f (.T3) - (rc - X2 -X3 + 1)0^(^3)) + 02(2:2) 

X((rc - 2:2 - 2:3 + l)0f (.X-3) + [Tc -X2 + XZ- 1)03 (Xs)))] - ^na(td) 

X [{n + rc + X2 - 1)02(X2)03 (xs) + T^r^ (02(x2)((4rfe + rc + 2:2 - 2:3 

-1)0^(0:3) - (rc + 2:2 + 2:3 - 1)0^(2:3)) + 0*2(2:2) ((rc + 2:2 + 2:3-1) 

x0f (2:3) - (rc + 2:2 -2:3 - 1)03(2:3)))] /i^«(2:2, 2:3, 61, 62)} , (31) 

The total decay amplitudes of the fifteen Be PV, VP decays can therefore be written 
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7r+ rP- 7r+ ~1 



-[/b.^;K«i+m:Kci]} , 

+ M^li'^^^C,] } cos , 

+[/Bci^;r''^«i + K'^^'^'^^i] } sin , 



A{B, ^ p+TfO) 



^(5, ^ p+rj) 



A{B, ^ p+r?') 



A{B, ^ /r^+r;) = V:,Vus fB. 



+ 

A{B^ ^ X*+?7') = K 



na 



cos0-M;^^^*"sin 



F;^*"''^sin</. + F;^^*%os0lai 



+ [m^:"-^" sin + M„^f *^ cos 0] } , 
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C. Bc^VV decays 



There are three kinds of polarizations of a vector meson, namely, longitudinal (L), 
normal (N), and transverse (T). The amphtudes for a Be meson decay to two vector 
mesons are also characterized by the polarization states of these vector mesons. The 
decay amplitudes ^A^'^^ in terms of helicities, for Be V{P2, £2)^(^3, Cg) decays, can be 
generally described by 

M^''^ - e;,{a)elMU9"' + W + i e'^'"''^ P2aP3p] , 

= ml Ml + m\M^e\{a = T) ■ e;{a = T) 

+iMTe''P^^eUa)elp{a)P2,P^p , (47) 

where the superscript a denotes the helicity states of the two vector mesons with L{T) 
standing for the longitudinal (transverse) component. And the definitions of the ampli- 
tudes Mi {i = L, N, T) in terms of the Lorentz-invariant amplitudes a, b and c are 

m% Ml = a e;{L) ■ el{L) + ^ el{L) ■ P, el{L) ■ P2 , 
Mn = a , 

ml Mt ^ ■ (48) 

We therefore will evaluate the helicity amplitudes Ml, Mn, Mt based on the pQCD 
factorization approach, respectively. 

For every component of the polarization, the corresponding Feynman amplitude can 
be written as the following form. 



/■I poo 

F^^ — SnCFinl^ / dx2dxz I b2db2bsdbs 
" Jo Jo 



X { [X2MX2)MX3) - 2r2r3 {{X2 + 1)02(2^2) + {X2 - 1)(I)1{X2)) 
X(pl{Xs)]Efa{ta)hfa{l - ^s, Xs, 63, 62) + E fa{tb)hfa{x2, 1 " 2:3, 62, ^3) 

X [(X3 - l)<p2{x2)Hxz) - 2r2r^<Pl{x2) ((xs - 2)0^(2:3) - x^cp\{x^))] ] , (49) 



^na = -^—^CfITiI^ J dX2 dx^ J 6iCi6i620?&2 



X {Ena{tc) [{Vc - X3 + 1)02(2^2)03(3^3) " '^2^3 (02(^:^2) ((Sr^ + X2 - X3 + 1) 

X03(a:3) - (rc -X2-X3 + 1)03(^^3)) + 0*2(a:;2)((rc - ^2 - X3 + 1)03(2^3) 

+ (rc -X2+X3- l)0^(a;3)))] /i^„(x2, X3, 61, 62) - Ka{^2: X3, 61, b2)Ena{td) 

X [{n + r,. + X2 - l)02(a;2)03(.«3) - r2r3 (02(.T2)((4r;, + rc + X2- X3 - 1) 

X03(a;3) - {tc + X2 + X3- 1)03 (^3)) + 0*2(x2)((re + X2 + X3- 1)03(^3) 
-(r, + X2-X3-l)0*3(x3)))]} , (50) 
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= SnCprn'^Bc / dx2dx^ / h2dh2h^dh^r2r^ 
" Jo Jo 



'0 ~'0 

X {hfa{l - X^,X2MM)Efa(ta) [{X2 + 1) (0^(a;2)03 (Xg) + 02 (a;2)03 (^^s) ) 
+ (X2 - 1)(0^(X2)0^(X3) + 02(3^2)03(^3))] + Ef^{h)hf,{x2, 1 - Xs, 62, ^s) 

X [(X3 - 2)(0^(a:2)0^(x3) + 0^(x2)0;^(x3)) - X3 (0^(a;2)0;^(x3) + 0^ (0:2)0? (0:3 ))i^l,) 



= TcCFrn^B / dx2 dx^ / bidbib2db 

3 'Jo Jo 



2?^2r3 



X {re [02(a;2)03(a;3) + 02(x2)03(a;3)] Ena{tc)hl^{x2, X^, 61, 62) 

-rb[0^(a;2)0?(a;3) + 0^(x2)0^(x3)]i?.a(td)/ii(x2,a;3,6i,62)} , (52) 

T 2 r /"°° 

= IQrcCFm^^ / dX2dx3 / 62C?&2&3C?&3 '^2'^3 

Jo Jo 

X - X3, X2, 63, b2)Efa{ta) [(^2 + 1) (0^(X2)03(X3) + 0^ (X2)0;^ (2:3) ) 

+ {X2 - 1)(0^(X2)0;^(X3) + 0^(X2)0^(X3))] + hfa{x2. 1 " 0:3, 62, b^)Efa{h) 

X [(X3 - 2)(0^(X2)0;;(X3) + 0^(X2)0?(X3)) - X3 (0^(x2)0?(x3) + 0^(x2)0^(x3))i^3,) 

^rfa = ^ — 7rC^m|^ / dx2dx3 / bidbib2db2r2n 
-J Jo Jo 

X {r, [0^(X2)03(X3) + 0^(x2)03(a;3)] Ena{tc)K^{x2, X3, 61, 62) 

-r,[0^(a;2)0^(x3) + 0^(x2)0?(x3)]F„a(td)/^i(a:2,X3,6i,&2)} • (54) 

For seven VV decays, considering all the polarization(if = L, N, T) contri- 

butions and the Feynman decay amplitudes as shown in Eqs. fl49H5^ . the total decay 
amplitude of these channels can be obtained directly, 

A1^(i?.-pV) = V:,V^,{[fsFfJ^a, + MCS^C,] 

-[fsFfitna, + M^^'^C,]] = , (55) 

+ [fBF;^f<^i + M2$"]C,] , (56) 
A^^(S. rV*+) = {/i..i^5>*"ai + m5>*"Ci} , (57) 



kaHhz, . T^*0 „+\ _ ./7^\AH/r, , „0n 



(58) 
(59) 

M'^iB, K*^p+) = v^M^'iB, ^ /r+p°) , (60) 



(61) 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, we will calculate the branching ratios ( and polarization fractions, rela- 
tive phases) for those considered thirty Be M2M3 decay modes. The input parameters 
and the wave functions to be used are given in Appendix |Al In numerical calculations, 
central values of input parameters will be used implicitly unless otherwise stated. 

TABLE I: The pQCD predictions of branching ratios(i? ii's) for Be PP modes. The dom- 
inant errors are induced from charm quark mass mc = 1.5 it 0.15 GeV, combined Gegenbauer 
moments of related meson distribution amplitudes(See Appendix [Al exphcitlv) . and the chiral 
enhancement factors mg = 1.4 it 0.3 GeV and ttiq = 1.6 it 0.1 GeV, respectively. 



Decay Modes 








Decay Modes 








(AS = 0) 


BR's{W-^) 






{AS 


= 1) 


BR's{W~^ 


) 




Be ^ TT+vrO 









Be- 








)^o:3("^o) 


Be TT+T] 


22.8tliime] 






Be- 


> K+r] 








Be ^ TT+T]' 


15.3^15 k: 






Be- 


> K+r]' 








Be ^ K+K° 


24.0l2;4(^^; 






Be- 




2.0l[Jj(me: 


-0.8 ("i 





For Bf. PP, PV, VP decays, the decay rate can be written as 

r = ^|^|^(fie ^ MsMs)^ (62) 

where the corresponding decay amphtudes A have been given explicitly in Eqs. (12011271) 
and Eqs. (13211461) . Using the decay amplitudes obtained in last section, it is straightforward 
to calculate the branching ratios with uncertainties as presented in Tables (lll lllll) . 



TABLE II: Same as Table U but for B^. PV modes. 



Decay Modes 








Decay Modes 






{AS = 0) 


BR's{W-'^ 


) 




{AS 


= 1) 


BR's{W-^ 


) 


Be - TT+^oO 


l-7To{me] 






Be- 




3.lt[Ji(me) 


^L5(«0lSi("^o) 


Be ^ K°K*+ 


1.8l°:I(me) 






Be- 




6.ltli{me) 


i2:9(«0tSJ(™o) 


Be TT+LO 


5.8tli{mc] 






Be- 




2.3iJi(me) 


t\i{ai)±0.l{mo) 



For Be VV decays, the decay rate can be written explicitly as, 

where |Pc| = |P2z| = iPszl is the momentum of either of the outgoing vector mesons. 
Based on the helicity amplitudes fllSl) . we can define the transversity amplitudes, 

Al = -^mlML, Aii=^V2mlMN, A^ = ^m%^/2{r^ - 1)Mt ■ (64) 
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TABLE III: Same as Table U but for VP modes. 



Decay Modes 
{AS = 0) 


BR's{W-'^) 


Decay Modes 
(A5 = 1) 


S/Z's(10-8) 


Be ^ P+TtO 

Be p^r] 

Be ^ P+V' 

Be ^ K*°K+ 


5.4l2-i(„i^)+0-9(«.)±0.0(mo) 


Be ^ K*07r+ 
Be ^ K*+7r0 
Be ^ K*+r? 
Be ^ K*+r/' 
Be ^ </.K+ 


i.6i°:fK)iO:?(a,)l°:S(mo) 

0.9+°:i(me)l°:^(ai)±0.0(mo) 
3.8±l.l(me)lJ-°(ai)±0.0(mo) 



for the longitudinal, parallel, and perpendicular polarizations, respectively, with the nor- 



malization factor C, = Y G'|,Pc/(167rm^^r) and the ratio r = P2 ■ P-i/{mM2 ' ^Ms)- These 
amplitudes satisfy the relation, 

i^Lr+iAi'+i^^i' = i (65) 

following the summation in Eq. ( l63|l . 

Since the transverse-helicity contributions manifest themselves in polarization observ- 
ables, we therefore define two kinds of polarization observables, i.e., polarization fractions 
{II, f\\, f±) and relative phases (0||, (f)±) as [53], 

r |^L{||,±)P , _ ^||(±) 

It should be noted that the final results of relative phases will plus one value, i.e., tt, due 
to an additional minus sign in the definition of Al- 

We also define another two quantities reflecting the effects of CP-violating asymmetries 



indirectly [53|, |54 



A^„ = ^ , A^, = , (67) 

where 0|| and 0_l are the CP-conjugated relative phases corresponding to (f)\\ and (f)±, 
respectively. 

With the complete decay amplitudes, by employing Eq. fl63|) and the input parameters 
and wave functions as given in Appendix \^ we will present the pQCD predictions for 
CP-averaged branching ratios, longitudinal polarization fractions and relative phases of 
the considered decays with errors as shown in Tables IIVI and |Vl 

Based on the pQCD predictions as given in Tables I -V, we have the following remarks: 

• Among considered pure annihilation Be PV/VP, VV decays, the pQCD predic- 
tions for the CP-averaged branching ratios for those AS* = processes are much 
larger than those of AS = 1 channels ( one of the two final state mesons is the 
meson ), which are mainly due to the large CKM factor |Kid/K(sP ~ 19. For 
Be 7r+7r°, p-p° decays, the contributions from uu and dd components cancel each 
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TABLE IV: The pQCD predictions of branching ratios{BR' s) and longitudinal polarization 
fractions(LPi<"s) for Bf, VV modes. 



Decay Modes 


BR'silQ-'^) 


LPF's (%) 


Be ^ 

Be p'^uj 
Be ^ K*°K*+ 




10.6i3:2(me)^2;i(a,) 
10.0tnK)l|^(a.) 


92.0l°i(me)l?i(a,) 


Be ^ K*°p^ 
Be ^ K*+/)0 
Be ^ K*+W 
Be ^ 0K*+ 


0.3t°:[](me)l°:°(a.) 


94.9j:[J-2(me)t?:°(a,,) 
94.9l°;2(me)l}i(a.) 
94.8l°;3(me)l}i(a.) 
86.4inK)l^;^(a.) 



TABLE V: The pQCD predictions of relative phases for Be — > VV modes. 



Decay Modes 


(All (rad) 


(/)_L (rad) 






Be ^ 

Be p'^UJ 

Be ^ K*°K*+ 


3.86i°:i(me)t°:?^(a.) 


4.43l°:}6,K)t°:?^(a.) 
3.76t°:J6(m,)tO:^8(a,) 






-7r/2 
-7r/2 


Be ^ K*V^ 

Be ^ K*+pO 

Be ^ K*+LJ 
Be ^ <^K*+ 


4.1lt°io'(m,)t°iO(a.) 
4.15t°ii(me)t°:2^(a.) 
3.80i°|^(m,)+°:^^(a.) 


4.20t°:J^(m,)t°:^?(a.) 

4.20t°:J^(m,)t°i?(a.) 

4.23l°:Ji(m,)tO:26(a,) 
3.89^°:?^(,n,);o.43(„^) 








-7r/2 
-7r/2 
-7r/2 
-7r/2 



other exactly and result in the zero branching ratios. In fact, these two channels 
are forbidden, even if with final state interactions. Simply, two pions can not form 
an s wave isospin 1 state, because of Bose-Einstein statics. Any other nonzero data 
for these two channels may indicate the effects of exotic new physics. 

There is no CP violation for all these decays within the standard model, since there 
is only one kind of tree operators involved in the decay amplitude of all considered 
Be decays, which can be seen from Eq. f[T2]) . 



The pQCD predictions for the branching ratios of considered B^. decays vary in 
the range of 10~ ( for Be ^ K K~^, K K*^ and p'^uo decays) to 10"^ ( for most 
AS" = 1 Be decays). The decays with the branching ratio of 10~^ can be measured 
at the LHC experiment Sq . 



As mentioned in the introduction, the authors of Ref. [38| studied many pure an- 
nihilation Be decays by employing the SU(3) flavor symmetry and the OGE model 
respectively, and presented their numerical estimates for the branching ratios of 
Be (f)K^,K K^,K and K K*^ decays. As a comparison, we show in 
Table VI the pQCD predictions and the results as given in Ref. [38] for relevant 
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TABLE VI: The pQCD predictions of branching ratios for and K^*' K^*^+ 

modes. As a comparison, the numerical results as given in Ref. |38|] are also listed in the last 
two columns. 



Channels 


pQCD Predictions 


SU(3) Symmetry 


OGE model 


Br(Bc ^ 0K+) 
Br(Bc ^K°K+) 
Br(Bc ^ K°K*+) 
Br(Bc ^ K*°K+) 
Br(Bc ^ K*°K*+) 


^■Qt],l{m,)t];l{a^) X 10-8 
2At',l{m,m{a.i) x 10"^ 
1.8+°:I(m,)l^:i(a,) X 10-7 
1.0±0.1(mc)lo:3(a») X 10"*^ 
l-0^[!:J(m.)i°|(a.) X 10-6 


0(10-7 ^ 10-8) 
0(10-6) 

0(10-6) 
0(10-6) 


5 X 10-9 
6.3 X 10-8 

9.0 X 10-8 

9.1 X 10-8 



channels. From Table VI, one can see easily that the pQCD predictions basically 
agree with the results obtained based on the SU (3) flavor symmetry. 

For Be {n'^ , p^){ri,ri') decays, the relevant flnal state mesons contain the same 
component uu + dd, they therefore have the similar branching ratios. The small 
differences among their branching ratios mainly come from the different mixing 
coefficients, i.e., cos0 and sin0. 

For Be K^T]'^'^ decays, however, one finds that Br(Bc K^t]') ~ 10 x Br{Bc 
K^rf), which is rather different from the pattern of Br{Bc tt^i]) ~ Br{Bc 
TT'^f]') and Br{Bc P^v) ~ Br {Be P^v')- This large difference can be un- 
derstood as follows: For the AS* = 1 processes, both rjg and rjg will contribute to 
Be — >■ K~^r] and K~^r]' decays but with an opposite sign for rjg and rjs term, as well 
as different coefficients. Which results in a destructive interference between rjq and 
ris component for B^ K'^rj, but a constructive interference for B^. — >■ K^i]' . This 



situation is very similar with that for the B —>■ Krj and Kt]' decays [55|, ISy, l57 



Unlike B^ K+r]^"^ decays, Br{Be K*+r]') ^ 45r(5, K*+r]) ~ 3.8 x 10-8. 
The reason is that both of them are mainly determined by the factorizable contri- 
butions of rjs term. 

For Be VV decays, we can find that (a) the branching ratios are in order of 
0(10-8 ~ 10-7) except for Br(Bc ^ K*°K*+) and Br(Be p+cj) ~ 10-6; (b) the 
longitudinal polarization fractions are around 95% within the theoretical errors ex- 
cept for Be — > (j)K*^ ( ~ 86%) and play the dominant role. 

According to the discussions in Ref. [3^, there are some simple relations among 
some decay channels in the limit of exact S'f/(3) fiavor symmetry. For Be —>■ PP 
decays, such relations are 

A{Be K%+) = V2A{Be K+n^) = XA{Be K+K^) , (68) 

where A = Ks/Kd ~ 0.2. For Be VP/PV and Be VV decays, the relations 



15 



read 



A{B^ K*%+) = V2A{B, K*+7c^) = XA{B, K*^K+), (69) 
A{B, p+K^) = V2A{B, p^K+) = XA{B, K*+K^), (70) 
{-IYA{B+ p+K*"") = {-lYV2A{B+ p'^K*+) = \A{B, K*+K^) (71) 

where £ = 0, 1, 2 ^. From our pQCD calculations, we notice that the first equality of 
each of the above relations (I68H71I) are valid in isospin symmetry. They hold exactly 
in our numerical calculations. The second equality of each relations are only valid at 
exact SU(3) symmetry thus they are violated at the order of SU(3) breaking effect 
in our calculations. 

• Since the LHC experiment can measure the B^ decays with a branching ratio at 10^^ 
level, our pQCD predictions for the branching ratios oi Be ^ K K^, K K*^ and 
p^uj decays could be tested in the forthcoming LHC experiments. 

• For most considered pure annihilation B^ decays, it is hard to observe them even in 
LHC due to their tiny decay rate. Their observation at LHC, however, would mean 
a large non-perturbative contribution or a signal for new physics beyond the SM. 

• It is worth of stressing that the theoretical predictions in the pQCD approach still 
have large theoretical errors induced by the still large uncertainties of many input 
parameters. Any progress in reducing the error of input parameters, such as the 
Gegenbauer moments and the charm quark mass m^, will help us to improve the 
precision of the pQCD predictions. 



V. SUMMARY 

In short, we studied the two-body charmless hadronic Be PP, PV/VP, VV decays 
by employing the pQCD factorization approach based on the kx factorization theorem. 
These considered decay channels can occur only via the annihilation diagram and they will 
provide an important testing ground for the magnitude of the annihilation contribution. 

The pQCD predictions for CP-averaged branching ratios, longitudinal polarization 
fractions and relative phases are displayed in Tables (i mv|) . From our numerical evaluations 
and phenomenological analysis, we found the following results: 

• The pQCD predictions for the branching ratios vary in the range of 10~^ to 10~^, 
basically agree with the predictions obtained by using the exact SU(3) flavor sym- 
metry. The Be — *■ IC^ and other decays with a decay rate at 10"^ or larger 
could be measured at the LHC experiment. 

• For Be PV/VP, VV decays, the branching ratios of AS* = processes are basi- 
cally larger than those of AS = 1 ones. Such differences are mainly induced by the 
CKM factors involved: Vud ~ 1 for the former decays while Vus ~ 0.22 for the latter 
ones. 



^ Here, since the longitudinal contributions dominate the Be K*^p^ decay, we use its longitudinal 
part (i.e., £ = ) to compare with the decay amplitude of Be K*'^K^ decay. 
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• Analogous to -B — > Kr]^'^ decays, we find Br (Be K^r]') ~ 10 x Br (Be K^r]). 
This large difference can be understood by the destructive and constructive inter- 
ference between the rjg and rjs contribution to the B^. —>■ K^r] and B^ —>■ K^r]' 
decay. 

• For Be VV decays, the longitudinal polarization fractions are around 95% except 
for Be (f)K*^ ( fi ~ 86%) and play the dominant role. 

• Because only tree operators are involved, the CP-violating asymmetries for these 
considered Be decays are absent naturally. 

• The pQCD predictions still have large theoretical uncertainties, induced by the 
uncertainties of input parameters. 

• We here calculated the branching ratios and other physical observables of the pure 
annihilation Be decays by employing the pQCD approach. We do not consider the 
possible long-distance (LD) contributions, such as the re-scattering effects, although 
they may be large and affect the theoretical predictions. It is beyond the scope of 
this work. 
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APPENDIX A: INPUT PARAMETERS AND DISTRIBUTION AMPLI- 
TUDES 

The masses (GeV), decay constants (GeV), QCD scale (GeV) and B meson lifetime 

are 



(/=4) 
'MS 


= 0.250, 


mw 


= 80.41, 


ms, = 


6.286, 




= 0.489, 




= 1.02, 


fct> 


= 0.231, 


/J = 


0.200, 


rriK* 


= 0.892, 


Ik* 


= 0.217, 


fx* 


= 0.185, 


nip = 


0.770, 




= 0.209, 




= 0.165, 




= 0.782, 




0.195, 


f 


= 0.145, 


< 


= 1.4, 


< 


= 1.6, 


m[J' = 


1.08, 


ml' 


= 1.92, 


rrib 


= 4.8, 




= 0.131, 


fK = 


0.16, 




= 0.46 ps 



For the CKM matrix elements, here we adopt the Wolfenstein parametrization for the 
CKM matrix, and take A = 0.814 and A = 0.2257, p = 0.135 and r/ = 0.349 @. 
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The twist-2 pseudoscalar meson distribution amp litude 4>p{P = vr, K) , and the twist-3 



ones (f)p and 0p have been parametrized as [58 



60 



ff 



(f)j,{x) = ^/L_ 6x(l - x) l+a^Cl^\2x-l) + a^C^'\2x-l) + a^C2'\2x~l)[A2) 
fp 



P^3/2, 



P^3/2, 



0p(x) 



2y/2K 



ff 



1 + ( 30r^3 - ^pI ] C]J\2x - 1) 



9 



1 + 6 ( 5?73 - ^r/gcus - - ^PpO^ ) (1 - lOx + lOx 



2v^ 



;i - 2x) 



(A3) 
(A4) 



with the Gegenbauer moments a\ = 0, af" = 0.17±0.17, = 0.115±0.115, af = —0.015, 
the mass ratio pTr{K) = '^■n{K)l'^Q^^ and Pr),(^) = 2mq(s)/mgg(ss), and the Gegenbauer 
polynomials C^(t), 



^(3t2-l) = l(3-30t2 + 35t^) 

2 8 



= 3t, (t) = ^(5t^-l) , (t) = ^(l-14t^ + 21t^) . (A5) 

In the above distribution amplitudes for kaon, the momentum fraction x is carried by the 



s quark. For both the pion and kaon, we choose 173 = 0.015 and u;3 = —3 [58|, |59| . 

The twist-2 distribution amplitudes for the longitudinally and tranversely polarized 
vector meson can be parameterized as: 



(^y{x) 



fy{x) 



3/5/ 



x(l — x) 



1 + 4yCl^\2x - 1) + a^^yC'^^'{2x - 1 



II r^3/2. 



x(l -x) 1 + afyCl^\2x - 1) + afyC^'\2x - 1) 



3/2/ 



(A6) 
(A7) 



Here fv and fy are the decay constants of the vector meson with longitudinal and 
tranverse polarization, respectively. The Gegenbauer moments have been studied ex- 
tensively in the literatures [giI, 62 |, here we adopt the following values from the recent 



updates [63|, l6J, |65 



a'l^. = 0.03 ± 0.02, a^l = a^^ = 0.15 ± 0.07, a^^, = 0.11 ± 0.09, a^^ = 0.18 ± 0.0^A8) 
aij^. = 0.04 ± 0.03, = a^^ = 0.14 ± 0.06, a^^, = 0.10 ± 0.08, a^^ = 0.14 ± 0.0fA9) 



The asymptotic forms of the twist-3 distribution amplitudes (py^ and are [42j: 



^^^-px - l)\ (t>v{x) = -^{2x - 1) 



2^6 
3/v- 
8^6 



(l + (2x-l)^), 



2^6 
3/y 
4^6 



(2x- 1). 



(AlO) 
(All) 
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APPENDIX B: RELATED HARD FUNCTIONS 



In this section, we group the functions which appear in the factorization formulae. 

The functions h in the decay amphtudes consist of two parts: one is the jet function 
St{xi) derived by the threshold re-summation [3], the other is the propagator of virtual 
quark and gluon. They are defined by 



hfa{x3,X2,b3,b2 



,171 , 



+d{h2 - h)H^^\^2MBM)U^2MBM) 



H^i\V^^MB^hm) 



1 - h)Hl,'\y/x2il - Xs)MbM)MVx2{1 - X3)MBb2 
+e{h2 - h)H'i\^X2{l - Xs)MB^h)MV X2{1 - X,)MbM) 

(1) 



X 



ITT TT\ 



M./f^)MbM 



> 



where 



(rc - X2)(l - X3) + r, 



2 

c ' 



Fd = rl - {1 - Vc - X2)x3 



(B2) 



(B3) 



SindH^^\z)=Jo{z) + iYo{z). 
The hard scales are chosen as 



ta = max{v^MB^, 1/62, 1/&3}, 
tb = max{Vl - xsMb,, I/&2, 1/&3}, 

= max{v/x2(l - X3)Mb,, ^/\{r, - X2)(l - X3) + r^lM^,, l/h, I/62}, 



td = max{^X2(l -X3)Mb, 



;i -r,-X2)x3|MB„ 1/61, 1/62}. 



(B4) 
(B5) 
(B6) 

(B7) 



They are given as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 

The St re-sums the threshold logarithms In^ x appearing in the hard kernels to all 
orders and it has been parameterized as 



2i+2T(3/2 + c) 

St{x) = ^ [x{l - x)Y 

V7rr(l + c) 



(B8) 



with c = 0.4±0.1. In the nonfactorizable contributions, St{x) gives a very small numerical 
effect to the amplitude [ggJ. Therefore, we drop St{x) in hna- 

The evolution factors Eja and Ena entering in the expressions for the matrix elements 
(see section III) are given by 



Efa{t) = a,{t)exp[-S2{t)-S3{t)], 

Ena{t) = as{t)exp[-SBit) - S'2(t) - 5*3 (t)] 162=63, 



(B9) 
(BIO) 
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in which the Sudakov exponents are defined as 



Ssit) 

S2{t) 



s r 



Ms \ 5 /-* dfi 



V2 



^^^^ h 



l/b2 



(Bll) 



iMsm, (B12) 



with the quark anomalous dimension 7^ = —as/n- Replacing the kinematic variables of 
M2 to M3 in 5*2 , we can get the expression for S^. The explicit forms for the function 
s{Q,b) are defined in the Appendix A in Ref. 3^ . 
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